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Abstract 

It is well known that the successful operation of cognitive radio (CR) between CR transmitter and CR 
receiver (CR link) relies on reliable spectrum sensing. To network CRs requires more information from spectrum 
sensing beyond traditional techniques, executing at CR transmitter and further information regarding the spectrum 
availability at CR receiver. Redefining the spectrum sensing along with statistical inference suitable for cognitive 
radio networks (CRN), we mathematically derive conditions to allow CR transmitter forwarding packets to CR 
receiver under guaranteed outage probability, and prove that the correlation of localized spectrum availability 
between a cooperative node and CR receiver determines effectiveness of the cooperative scheme. Applying our 
novel mathematical model to potential hidden terminals in CRN, we illustrate that the allowable transmission region 
\Q • of a CR, defined as neighborhood, is no longer circular shape even in a pure path loss channel model. This results 

in asymmetric CR links to make bidirectional links generally inappropriate in CRN, though this challenge can be 
alleviated with the aid of cooperative sensing. Therefore, spectrum sensing capability determines CRN topology. For 
multiple cooperative nodes, to fully utilize spectrum availability, the selection methodology of cooperative nodes 
is developed due to limited overhead of information exchange. Defining reliability as information of spectrum 
availability at CR receiver provided by a cooperative node and by applying neighborhood area, we can compare 
sensing capability of cooperative nodes from both link and network perspectives. In addition, due to dynamic 
network topology lack of centralized coordination in CRN, CRs can only acquire local and partial information in 
limited sensing duration, robust spectrum sensing is therefore proposed to ensure successful CRN operation. Limits 
of cooperative schemes and their impacts on network operation are also derived. 
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I. INTRODUCTION 

OGNITIVE radios (CR) JUS, having capable of sensing spectrum availability, is considered as a 
^ ] promising technique to alleviate spectrum scarcity due to current static spectrum allotment policy . 

c3 ■ Traditional CR link availability is solely determined by the spectrum sensing conducted at the transmitter 
(i.e. CR-Tx). If the CR-Tx with packets to relay senses the selected channel to be available, it precedes 
this opportunistic transmission. To facilitate the spectrum sensing, at time instant t n , we usually use a 
hypothesis testing as follows. 

#1 : Y = I + N 

(1) 

H :Y = S + I + N 

where Y means the observation at CR-Tx; S represents signal from primary system (PS); I is the 
interference from co-existing multi-radio wireless networks; N is additive white Gaussian noise (AWGN). 
They are all random variables at time t n . We can conduct this hypothesis testing in several ways based 
on different criterions and different assumptions 
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1) Energy Detection ||9l- [fT3l : Energy detection is widely considered due to its simple complexity and 
no need of a priori knowledge of PS. However, due to noise and interference power uncertainty, 
the performance of energy detection severely degrades [f8l [|20l , and the detector fails to differentiate 
PS from the interference. 

2) Cyclostationary Detection lfT6l - lfT8ll : Stationary PS signal can be exploited to achieve a better and 
more robust detector. Stationary observation coupled with the periodicity of carriers, pulse trains, 
repeating spreading code, etc., results in a cyclostationary model for the signal, which can be 
exploited in frequency domain by the spectral correlation function lfl4"llfl"5l . with high computational 
complexity and long observation duration. 

3) Locking Detection lT71lfP9l : In practical communication systems, pilots and preambles are usually 
periodically transmitted to facilitate synchronization and channel estimation, etc. These known 
signals can be used for locking detection to distinguish PS from noise and the interference. However, 
locking detection requires more a priori information about PS, including frame structure, modulation 
types, and coding schemes, etc. 

4) Covariance-based Detection ||20l - [|24l : Because of the dispersive channels, the utility of multiple 
antennas, or even over-sampling, the signals from PS are correlated and can be utilized to differ- 
entiate PS from white noise. The existence of PS can be determined by evaluating the structure of 
covariance matrix of the received signals. The detector can be implemented blindly [24] by singular 
value decomposition, that is, it requires no a priori knowledge of PS and noise, but needs good 
computational complexity. 

5) Wavelet-based Detection ||25l : The detection is implemented by modeling the entire wideband signals 
as consecutive frequency subbands where the power spectral characteristic is smooth within each 
subband but exhibits an abrupt change between neighboring subbands, at the price of high sampling 
rate due to wideband signals. 

However, above spectrum sensing mechanisms, focusing on physical layer detection or estimation at 
CR-Tx, ignore the spectrum availability at CR receiver. We could illustrate the insufficiency of traditional 
spectrum sensing model, especially to network CRs. Due to existence of fading channels and noise 
uncertainty along with limited sensing duration [4], even when there is no detectable transmission of PS 
during this venerable period, the receiver of this opportunistic transmission (i.e. CR-Rx) may still suffer 
from collisions from simultaneous transmission(s), as Fig. Q] shows. The CR-Rx locates in the middle of 
CR-Tx and PS-Tx and PS activities are hidden to CR-Tx, which induces a challenge to spectrum sensing. 
We can either develop more powerful sensing techniques such as cooperative sensing [12611 - 11321 to alleviate 
hidden terminal problem, or a more realistic mathematical model for spectrum sensing what we are going 
to do hereafter. 

The organization of this paper is as follows. We elaborate a realistic definition of link availability and 
system model in Section II and present general spectrum sensing with/without cooperation in Section III 
an IV respectively. In Section V, realistic operation of CRN is suggested to investigate impacts of spectrum 
sensing and cooperative scheme on network operation. Numerical results and examples are illustrated in 
Section VI. Finally, conclusions are made in Section VII. 

II. System Model 

From a viewpoint of information theory, spectrum sensing can be modeled as a binary channel that 
transmit CR link availability (one bit information in link layer) to CR-Tx with transition probabilities 
representing spectrum sensing capability, probability of missing detection and probability of false alarm. 
Therefore, traditional spectrum sensing mechanisms could be explained by a mathematical structure of 
defining link availability. 
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Definition 1. CR link availability, between CR-Tx and CR-Rx, is specified by an indicator function 



1, CR link is available for opportunistic transmission 
0, otherwise 



Definition 2. CR-Tx senses the spectrum and determines link availability based on its observation as 

^ Tx I 1, CR link is available for transmission at CR-Tx 
1 0, otherwise 

Lemma 1. Traditional spectrum sensing for CR link suggests l hnk = l Tx . 

The definition of link availability is pretty much similar to the clear channel assessment (CCA) of 
medium access control (MAC) in the IEEE 802.11 wireless local area networks (WLAN), or the medium 
availability indicator in 11471 . We may have a correspondence between link availability in dynamic channel 
access of cognitive radio networks (CRN), and the CCA in MAC of WLAN. 

As we explain in Figure \T\ and/or take interference into testing scenario, we may note that Lemma 
Q] is not generally true. To generally model spectrum sensing, including hidden terminal scenarios, we 
have to reach two simultaneous conditions: (1) CR-Tx senses the link available to transmit (2) CR-Rx 
can successfully receive packets, which means no PS signal at CR-Rx side, nor significant interference 
to prohibit successful CR packet reception (i.e. beyond a target SINR). In other words, at CR-Rx, 

P(j Tx 

SINRcR-Rx = 75 — 5 — 5 > Vantage (2) 

" PS + "l + -*/V 

where r] outag e is the SINR threshold at CR-Rx for outage in reception over fading channels, Pcr-tx 
is the received power from opportunistic transmission from CR-Tx, Pp$ is the power contributed from 
PS simultaneous operation for general network topology such as ad hoc, Pj is the total interference 
power from other co-existing radio systems ll38l . and Pn is band-limited noise power, with assuming 
independence among PS, CR, interference systems, and noise. 

Based on this observation, CR link availability should be composed of localized spectrum availability 
at CR-Tx and CR-Rx, which may not be identical in general. The inconsistency of spectrum availability 
at CR-Tx and CR-Rx is rarely noted in current literatures. However, this factor not only suggests spatial 
behavior for CR-Tx and CR-Rx but also is critical to some networking performance such as throughput 
of CRN, etc. Il33l developed a brilliant two-switch model to capture distributed and dynamic spectrum 
availability. However, [|33l focused on capacity from information theory and it is hard to directly extend 
the model in studying network operation of CRN. Actually, two switching functions can be generalized as 
indicator functions to indicate the activities of PS based on the sensing by CR-Tx and CR-Rx respectively 
I1471 . Generalizing the concept of [133111341 to facilitate our study in spectrum sensing and further impacts 
on network operation, we represent the spectrum availability at CR-Rx by an another indicator function. 

Definition 3. The true availability for CR-Rx can be indicated by 

^ Rx I 1, CR link is available for reception at CR-Rx 
1 0, otherwise 

Please note that the activity of PS estimated at CR-Rx in 11331 may not be identical to l Rx . That is, 
even when CR-Rx senses that PS is active, CR-Rx may still successfully receive packets from CR-Tx if 
the received power from CR-Tx is strong enough to satisfy (O. We call this rate-distance nature ||36l that 
is extended from an overlay concept 041 . Therefore, we consider a more realistic mathematical model 
for CR link availability that can be represented as multiplication (i.e. AND operation) of the indicator 
functions of spectrum availability at CR-Tx and CR-Rx to satisfy two simultaneous conditions for CR 
link availability. 
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Proposition 1. l lmk = 

To obtain the spectrum availability at CR-Rx (i.e. 1^) and to eliminate hidden terminal problem, a 
handshake mechanism has been proposed by sending Request To Send (RTS) and Clear To Send (CTS) 
frames. However, the effectiveness of RTS/CTS degrades in general ad hoc networks 11351 . Furthermore, 
since CRs have lower priority in the co-existing primary/secondary communication model, CRs should 
cherish the venerable duration for transmission and reduce the overhead caused by information exchange 
and increases spectrum utilization accordingly. Therefore, the next challenge would be that l Rx cannot be 
known a priori at CR-Tx, due to no centralized coordination nor information exchange in advance among 
CRs when CR-Tx wants to transmit. As a result, general spectrum sensing turns out to be a composite 
hypothesis testing. In this paper, we introduce statistical inference that is seldom applied in traditional 
spectrum sensing to predict/estimate spectrum availability at CR-Rx and to regard it as performance lower 
bound in general spectrum sensing. 

Further examining Proposition [H we see that prediction of l Rx is necessary when l Tx = 1, which is 
equivalent to prediction of \ hnk . In this paper, we model 1^ when l Tx = 1 as a Bernoulli process with 
the probability of spectrum availability at CR-Rx Pi(l Rx = l\l Tx = 1) = a. The value of a exhibits 
spatial behavior of CR-Tx and CR-Rx and thus impacts of hidden terminal problem. If a is large, CR-Rx 
is expected to be close to CR-Tx and hidden terminal problem rarely occurs (and vise versa). 

III. General Spectrum Sensing 

The prediction of l Rx at CR-Tx can be modeled as a hypothesis testing, that is, detecting l Rx with a 
priori probability a but no observation. To design optimum detection, we consider minimum Bayesian 
risk criterion, where Bayesian risk is defined by 

R = wPr(l link = 0\l Tx = 1)P F + p Y (l link = l\l Tx = 1)P M (3) 

In ©, P F = p r (i /mfe = l\l link = 0,l Tx = 1), P M = p r (i tofc = 0\l link = l,l Tx = 1), and w > 
denotes the normalized weighting factor to evaluate costs of Pp and P M , where l hnk represents prediction 
of l lmk . We will show that the value of w relates to the outage probability of CR link in Section V. 

Since l Rx is unavailable at CR-Tx, we have to develop techniques to "obtain" some information of 
spectrum availability at CR-Rx. Inspired by the CRN tomography ||39lll40l , we may want to derive the 
statistical inference of l Rx based on earlier observation. It is reasonable to assume that CR-Tx can learn 
the status of l Rx at previous times when l Tx = 1, which is indexed by n. That is, at time n, CR-Tx 
can learn the value of l Rx [n — 1], l Rx [n — 2], . . .. In other words, we can statistically infer l^n] from 
l Rx [n — 1], l Rx \n — 2], . . ., l Rx [n — L], where L is the observation depth. This leads to a classical problem 
from Bayesian inference. 

Lemma 2. Through the Laplace formula 143], the estimated probability of spectrum availability at CR-Rx 
is 

- N + 1 

a = (4) 

L + 2 v 

where N = ^f =1 l Rx [n - I]. 

Proposition 2. Inference-based spectrum sensing at CR-Tx thus becomes 

{1 Tx if n > w 
' y — w+1 
0, otherwise 

where a is in (I?]). 

Proof: Since the optimum detector under Bayesian criterion is the likelihood ratio test ||42)| . we have 

w (6) 



Pr( 




= l\l Tx 


= 1) 


_ > 


- 1 Cq\ — Coo 


Pr( 


'■]_Rx 


- 0\l Tx 


= 1) 


1 Ot ^Unk 


=0 Cio — C\i 
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where CV, denotes the cost incurred by determining l hnk = j when l hnk = i. According to Bayesian risk 
in ©, we have Cn = C 00 = and C i/Cio = w. Rearranging the inequality, we obtain the proposition. 

■ 

Remark. CR-Tx believes CR link is available and forwards packets to CR-Rx if the probability of spectrum 
available at CR-Rx a is high enough. Otherwise, CR-Tx is prohibited from using the link even when CR-Tx 
feels free for transmission because it can generate unaffordable cost, that is, intolerable interference to 
PS or collisions at CR-Rx. 



IV. General Cooperative Spectrum Sensing 

A. Single Cooperative Node 

Spectrum sensing at cooperative node, which can be represented l Co , is to explore more information 
about l Rx and therefore alleviates hidden terminal problem. We can use Fig. [2] to depict the scenario. In 
case the existence of obstacles, 1^ is totally orthogonal to l Tx . l Co is useful simply because of more 
correlation between l Co and 1^. From above observation, we only care about correlation of l Rx and 
l Co when l Tx = 1 and assume 

Pr(l Co = l|l fe = l,l T:E = l)=/? 
Pr(l Co = 0|l fe = 0,l T:E = l)=7 

Thus the correlation between l Co and l Rx , p, and corresponding properties become 



>/a(l-a)(/? + 7-l) (7) 
^ +(!_«)(!_ 7 )) (a( i _ (3) + (i _ a)7 ) 

Lemma 3. p is a strictly concave function with respect to a G (0,1) if 1 < [3 + 7 < 2 but a strictly 
convex function if < /3 + 7 < 1. In addition, l Co and 1^ are independent if and only if p = 0, i.e, 

13 + 7=1 

Proof: Let 

x / a(l-a) 
/(«) = 



{a(3 + (1 - a)(l - 7 ))(a(l - P) + (1 - a)j) 
Then taking first order and second order differentiation with respect to a, we have 

f/(n) = -a 2 /3(l-/3) + (l-a) 2 7(l-7) 

1 1 ' 2(a(l - a)y/\al3 + (1 - a)(l - 7 )) 3 / 2 (a(l - 0) + (1 - a) 7 ) 3 / 2 

f"(a) = K(a)[-I0a 2 (l - a) 2 /3 7 (l - ^(1 - 7) ~ 4a 4 (l - a) 13(1 - (3){(37 + (1 - P)(l ~ l)) 
+ (3 - 4«)« 4 /? 2 (l - (3) 2 - 4a(l - a) 4 7 (l - 7)(/?7 + (1 - /3)(1 - 7)) 
+ (-l + 4a)(l-a) 4 7 2 (l- 7 ) 2 ] 

where if (a) > for a, (3,7 e (0, 1) and (3 + 7 ^ 1. In addition, 0y + (1 - /3)(1 - 7) - /3(1 - /3) = 
/?(/? + 7 - 1) + (1 - /3)(1 - 7) = py + (1 - - p - 7) > for p,y e (0, 1). Similarly, we have 
Py + (1 — /3) (1 — 7) > 7(1 — 7) for P, 7 G (0, 1). Therefore, combining the second and the third terms, 
and the fourth and the last terms in the bracket of f"{a), we have 

f(a) < K(a)[—10a 2 (l - a) 2 Py(l - p)(l - 7) - a 4 p 2 (l - p) 2 - (1 - «) 4 7 2 (1 - 7 ) 2 ] < 
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if a G (0, 1). Since p = {(3 + 7 — l)/(a), we prove the first statement of the lemma. For the second 
statement, obviously, p = if l Co and 1^ are independent. Reversely, if p = 0, i.e., /3 + 7 = 1, we have 

1 

Pr(l Co = l\l Tx = 1) = J2 Fr ( lRx = s \ lTx = l)Pr(l C ° = l\l Rx = s,l Tx = 1) 

= (1 - a)(l - 7) + a/5 = (1 - a)/3 + a(3 = (3 

= Pi(l Co = l\l Rx = l,l Tx = 1) 

= 1 - 7 = Pr(l Co = l\l Rx = 0, l Tx = 1) 

Similarly, we can show Pr(l Co = 0\l Tx = 1) = Pr(l Co = 0\l Rx = l,l Tx = 1) = Pr(l Co = 0\l Rx = 
0,l Tx = 1) and complete the proof. ■ 
By statistical inference, CR-Tx can learn statistical characteristic of l Rx and l Co , i.e., {a, /3, 7}, by 
previous observations. From a viewpoint of hypothesis testing, we would like to detect \ Rx with a priori 
probability a and one observation l Co , which is the detection result at the cooperative node. In addition, 
probability of detection and probability of false alarm at the cooperative node are (5 and 1 — 7 respectively. 

Proposition 3. Spectrum sensing with one cooperative node becomes 

l Tx , if a > max{ai, a^} 

Vink _ ) l Tx l Co , ifa 2 <a<a l7 p>0 

l Tx lC°^ if an < a < a 2 , p < 

0, if a < min{a;i, a%) 

where l Co is the complement of l Co , a% = wj/(l — /3 + wj) and a 2 = w(l — j)/(/3 + w(l — 7)). 

Proof: The likelihood ratio test based on observed signal l Co can be written as follows. For l Co = 1, 



For l Co = 0, 



Pr(l Co 




= 1, 


-jTcc 


1) 









l-a) 


Pr(l Co 




= 0, 


| 


1) 


1 




^1/ £Zinfc_Q 


a 


Pr(l Co 


= 0|l to 


= 1, 




1) 


1 






[l-a) 


Pr(l Co 


= 0\l Rx 


= 0, 




1) 




1 


< 


a 



Rearranging the above inequalities, we obtain the proposition. ■ 
It is interesting to note that cooperative spectrum sensing is not always helpful, that is, it does not always 
further decrease Bayesian risk. We see that if a is large (greater than max{«i, a 2 }), that is, hidden terminal 
problem rarely occurs because either CR-Rx is close to CR-Tx or CR-Tx adopt cooperative sensing to 
determine l Tx , prediction of 1^ is unnecessary at CR-Tx. On the other hand, if a is small (less than 
minjai, a 2 }), CR-Tx is prohibited from forwarding packets to CR-Rx even with the aid of cooperative 
sensing. 

In the following, we adopt minimum error probability criterion (i.e., w = 1) and give an insight into the 
condition that cooperative sensing is helpful. Although we set w = 1, we do not lose generality because 
we can scale a priori probability a to a/(a + w(l — a)) as w ^ 1. Applying Lemma |3] and the fact that 
p\a=a 1 c — p\a=a 2 c when w = 1, we can reach the following corollary. 

Corollary 1. If we adopt minimum error probability criterion, spectrum sensing with one cooperative 
node becomes 

Wnk _ J ifei/a] 1 ^ if \p\ < * 9) 

l [p>0] l Co + l [p <o]l Co )l T *, if\p\ >* 



7 



where lr s i is an indicator function, which is equal to 1 if the statement s is true else equal to 0, and 



^2(07 +(!-/?)(! -7)) 

Remark. The effectiveness of a cooperative node only depends on the correlation of spectrum availability 
at CR-Rx and the cooperative node. If the correlation is low, information provided by the cooperative 
node is irrelevant to the spectrum sensing which degenerates to ([5]). 

By establishing a simple indicator model in link layer, we mathematically demonstrate the limit of a 
cooperative node in general spectrum sensing. It is natural to ask what will happen for multiple cooperative 
nodes and how to compare sensing capability among cooperative nodes. In the following, we provide 
metrics to measure sensing capability of cooperative nodes from link and network perspectives. 



B. Preliminaries 



Before exploring multiple cooperative nodes, we introduce notations and properties to systematically 
construct relation between joint probability mass function (pmf) and marginal pmf of spectrum availability 
at cooperative nodes. We first define notations in the following. 

Definition 4. For an m x n matrix A and two n x 1 vectors u and v, 



A[i,*] : 
A T : 

pJlR . 
Omxn • 

u© v : 
||u|| p : 



the ith row of A 
A T [i,j)=A\j,i) 
A RR [i, *] = A[m - i + 1, ★] 

am = o 

u v[z] = u[z]v[i] 

own 17 * 



Let 



q(1) 

m,k 

p(0) 



[An An 



A[*, j] : 


the jth column of A 


A RC. 


A RC {*,j] = A[*,n-j' + l] 


lmxn • 


A[i,j] = l 


In: 


n x n identity matrix 


u ^ v : 


u[i] < v[i] 


l[u>0] : 


l[u>0][«] = l[u[t]>0] 


A k,k\ 


< n < k 




■■ (AD T f 



where A™ , is a ( k ) x ( k ) matrix, A 
and for 1 < n < k — 1 



o 

m,k 



-lx 



0<m<A:, A^ fc = 1; 



,0 <m < k-1, A 



k 

k,k 



(10) 

(11) 

(12) 

= 1, 



m,k 



Xl 



Afe k 



An 

^m,k~l 



An 

-^-m-l.fe-l 
a n-1 

^Nn~l,fc-1 



1 < m < k — 1 



s = 0,1. In the following, we show properties of 



The role of A^ (or A l m k , < m < k) is to specify arrangements of joint pmf and marginal pmf such 

' (s) 

that their relation can be easily established by G m k , 
Ka,k and G%. 
Lemma 4. Let I m ,k(j) 



{z|A^ fe [i, j] — 1, 1 < i < k}, I < j < ( ) and we have 

\Zm,k{j)\ = m 

Zm,kU) = Zm,k(l) if and only if j = I 
where \I\ denotes number of elements in the set I. 



(13) 
(14) 
(15) 
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Remark. Since A^ fc is a k x (^) matrix, from rf7?l ), rf7?l) . we conclude that A^ fc [*, j], 1 < j < ( k ) 
contains all possible permutations of m ones and k — m zeros. 

Lemma 5. Let S m = YhLo (*)• Let G m,k and Q_m,k be S m x S m matrices, G^ k and G^ k be (2 k -S m ) x S m 



matrices, and G 



m,k 



p (s) 



1 

Ofcxl 



,s = 0,l. 

llxfc 



0/ 



0/ 



xl 



xi 



1 ixfM 



L 




A 1 







'CO 


°(i)x* 

r>(o) _ 

^m,fc ~~ 











Ma' 

A 1 



An 



are nonsingular and their inverse matrices become 

1 — llxfc • • • 



Ofcxl 

O/fe 



h 



"l) 1 ^//.fc 



ix r 



xl 







I 



G) 



C)xl °( 



)xfc 







(-^Ma 

("l) 1+m A^ fc 

(_l)n+m A n^ 



frK 1 ) \ —RR 



Let p^Ji] = Pr(lf° = A^[l, i], . . . , if = A^[*, z]|l fc = s, 1 T * 
and 1 < i < ( ), lf° denotes spectrum availability at the ith cooperative node, and let 

T 



(16) 



(17) 



= A^ fc [fc,z]|l te = s, l Tx = 1), where s = 0, 1, < m < fc 



»(*) 



(p$) T (p5 



to' 



(s) NT 



Therefore, characterizes joint pmf of spectrum availability at k cooperative nodes, . . . , l k °. 



Co 



(18) 



I Co 



Similarly, we divide Pj^ into two parts. Let P^' fc and P^ ; fc be S m x 1 vectors, ~P)^ k and P 



(i) 



,(0) 



j(0) 

m,k 



be 



(2 fc — S* m ) x 1 vectors, and P k 

qSfcbl = Pr(i£° 



fp< a > NT /p 

v m,kJ \—m,k 



O, . . . , Xju 



0) \T 
i,k) 

s\l Rx = 



, s = 0, 1. In addition, let 

s, l Tx = 1, {A;i, . . . , k m } = T m ,k(j)) 



s = 0, 1, 1 < m < k and 1 < j < ( ), which specifies the mth order marginal pmf of 1^ 
Arrange them into a vector form, we define 



Co 



Please note that q« = Q^[2 fe ] = P)?[2*] and q$ = Q^[2 fc ] = P™[1]. 



Q (5) fc 
rofc 



i (q5) T (qg) T 



C W \T 



i T 



,s = 0,l 



(i)rofci 



»(°) rofci 



>(0)[ 



Lemma 6. Marginal and joint pmf of spectrum availability at k cooperative nodes satisfy G m k ±. k 



W P W 



or 



(1) r»(l) 



p(0) 

i-fc,_fsT 



/p(0) \-l/ n (0J pTW 



l(0) 



7t(°) t^(°) 



iOo 
> J-fc • 



(19) 



Q (s) fc 



(20) 
(21) 
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Lemma 7. Pj^ provides equivalent information to Q^' k , s — 0, 1, or specifically, 



where 



cW + Q«[2*]bW ia = 0,l 



rr (1) hin (1) 



,(0) 



:-iYb 



(22) 



C. Multiple Cooperative Nodes 

Assume there are K cooperative nodes with corresponding spectrum availability if °, if °, 
their joint pmf conditionally on l Rx = s and l Tx = 1, Pj£ , s = 0, 1. 

Proposition 4. Spectrum sensing with multiple cooperative nodes becomes 

2 K K 

V ink = lTx r[j] "Q [A ^, j]lf © (1 - A^j})! 

3=1 i=l 

and 



. , 1 C K ° and 



Col 
k J 



^minHl-^P^^aP^i]} 



(23) 



(24) 



i=l 



where T — 1. _m s^ro^m © denotes OR operation. 

\aP K — w(l— a)P K >0\ 

Proof: By the likelihood ratio test, we have 

P^[z] i iirafc =i«;(l - a) 



> 
< 



Therefore, the optimum detector becomes 



r[z] — i m |ico =A ^.|- lii ] i „. )1 oo =A ^[ JC( , 



1 [aP«[i]-«;(l- a )p(, 0) [i]>0] 



With the result and combining binary arithmetic we obtain (1231) . In terms of Bayesian risk, 



Ro P t = E M 1 - tt ) P * M Pr ( iK " fc = il 1 ^ = A Ml, 



l£° = A^[AT, i), l Hx = 0, 1 T:E = 1) + 



I Tie 



i=l 



aP^[i\ Vi(l lmk = 0|lf° = A^[l,z], • • • , 1° K ° = A^K, i], l Hx = 1, l Tx = 1)] 

: £ Hi - «)Px[W] + oP^Ma - r[i])] 

i=l 

^minMl-^P^aP^]} 



i=l 



We consider two cooperative nodes insightfully to understand how multiple cooperative nodes improve 
performance of spectrum sensing. For spectrum availability at two cooperative nodes if °, if °, let = 

[fli f32\~ and qfl = [71 72] T ■ In addition, their joint probability is specified as follows. When l Rx = 1, 
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PS is likely inactive and then lf° and if" are independent. On the other hand, when l Rx = 0, if and 



1 2 ° are correlated with correlation p 12 . With the constraints on q^ u 2 , we have 



(0) 



,(0) 



[ 7l72 + A (l- 7l ) 72 -A 7l (l- 72 )-A (l- 7l )(l- 72 ) + A] : 



(25) 



where A = a/ 7 i 72 (1 - 7 i)(l - 7 2)pi2- 

1) Independent (p l2 = 0): In case lf° and 1 2 ° are conditionally independent, it leads to conventional 
assumption in cooperative spectrum sensing. 

Proposition 5. For two cooperative nodes with independent spectrum availability, spectrum sensing 
becomes 



i link 



i=l 



^l[Pi>o]l?° + l[p l <o]l l Co ), 



10, 



2 

i=X 



lfc>o]lf° + l 



[ft<0]- 



\Co\ 



if a> a (4) 
if "(3) < a < ct( 4 ) 
if a (2) < a < a (3 ), k 
if < a < a (2) 
if a < a (i) 



arg maxj M R (i) 



(26) 



where 



V«(l),«(2),a(3),«(4)) 



„c „,c „,c 



a? = wPf> [i]/(P 2 1} [i] +^Pf [i]), M fl (i 



,(i)r 



(«4 > °€ 
(of, af 

«f 
(of, of 

= 1 



Q 



a 



c 



a 



c 



a 



«4 



a. 



a 



a 



a 



a: 



a: 



a 



if Px > 0,p 2 > 0,5f > 5t 
if Px >0,p 2 > 0,8f < <5 + 
if Px > 0,p 2 < 0,5+ > 5 2 
if px > 0,p 2 < 0,5 + < 5 2 
if px < 0, p 2 > 0, <5r > 6% 
if Px < 0,p 2 > 0,5f < 5+ 
z/pi < 0,p 2 < 0, <5f > &f 
|f px < 0, p 2 < 0, <ff < 5 2 



(27) 



[ ft >o]^ + l[ Pl <o]^ . <V = 7iA/((l - 7i)(l 
(1 — 7i)(l — A)/(7iA)' Pi denotes the correlation between lf° and l Rx . 

We list all valid orders of (af , a^, a^, a^) in (ITTT) and the the spectrum sensing is determined according 
to the value of a by arguing (|23l) . We observe that a is high (ate) < a < au)), any one of cooperative 
nodes helps the spectrum sensing, which leads the spectrum sensing to OR operation. However, when 
a is low (a(i) < a < «( 2 )), CR-Tx requires more evidence to claim available link and the spectrum 
sensing becomes AND operation. In addition, it is interesting to note that there exists a region of a 
(«( 2 ) < a < a ( 3 )) such that CR-Tx only depends on one of two cooperative nodes, which motivates us to 
define a metric or a measure to evaluate cooperative nodes. 

Definition 5. Reliability of a cooperative node is measured by M R , if ° is said to be more or equally 
reliable than (to) lf° if M R (i) > M R (j), which is denoted by lf° > if". 

Proposition 6. For K cooperative nodes with independent spectrum availability, without loss of generality, 
we assume pi > for i = 1, . . . , n and p. t < for i = n + 1, . . . , K. Spectrum sensing becomes 
l hnk = l Tx l\ s ], where 



£ logM fl (z) >log 

ieC+uC- 



a 



a 



' 1=1 x ' «=n+l x 



(28) 



and C + = = 1, % = 1, . . . , n}, C~ = {i\lf° = 0,i = n + 1, . . . , K} 
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Proof: Since 1\ 



Co 



1 Co 
L 2 > 



1^° are independent, the likelihood ratio test becomes 



-7* 



n 

i€{l,...,n}\C+ 



l-A 

7* 



n 



l-A 



Ti 



n 



A 



A7i 



A)(l"7i) 



n 



M, 



iec+uc- 



n 

ie{l,...,n} 



n 

i6{l,...,n} 

i-A 



A 



7i 



n 



7i 



n 

ie{n+l,.. 



(i-A)(i-7*) 

A7i 
A i"^=iu;(l 



1-7, 



n 

ie{n+l,.. 



A 



,K} 



li 



a 



a 



Taking logarithm at both sides and rearranging the formula, we complete the proof. ■ 
We observe that reliability Mr{i) is used to quantify the information of l^ 1 provided by the zth 
cooperative node. Please note that if lf° is independent of l Rx , M R {i) = 1 and lf° is irrelevant to the 
spectrum sensing. Therefore, reliability can imply sensing capability, that is, one cooperative node with 
higher reliability has better sensing capability. Reliability can thus serve a criterion to select cooperative 
nodes when number of cooperative nodes is limited due to appropriate overhead caused by information 
exchange. Specially, if there are K equally reliable cooperative nodes, each cooperative node provides 
equal amount of information about l Rx and the spectrum sensing rule turns out to be Counting rule. This 
is a generalization from identically and independently distributed (i.i.d.) observations 11371 in conventional 
distributed detection to equally reliable observations. 

2) Correlated (p\i ^ 0): When there exists correlation between spectrum availability at cooperative 
nodes, the joint probabilities Pr, are shifted by A, as in (|25l) . For example, if the correlation is positive 
ase while af and a% decrease. If the correlation increases, eventually, the order of 
will switch and the spectrum sensing in (T26l) will change accordingly. However, the 



c 



a\ and a 4 



a 



correlated case becomes tedious and we consider a simple but meaningful example, where cooperative 



nodes have symmetric error rates (i.e. A = id and reliability becomes l[ Pi >o]A 
with the aid of (1231) and (l25l) . the spectrum sensing can be easily derived. 



-\Pi< 



oj(l-A)- Similarly, 



Proposition 7. For two cooperative nodes with correlated spectrum availability and symmetric error rate 
satisfying p x > 0, p 2 > 0, and lf°[> lf°, the spectrum sensing would be ([26]) with modifications according 
to A. 

if a (2 ) <a< a (3) , A < (1 - 2A)A 
ifap) <a< «( 4 ), A < (1 - 2A)A 
i/a ( i) < a < a (2) , A > (2/3 2 - 1)(1 - A) 
ifa {2) < a < «(3), A > (2A - 1)(1 - A) 



I link 



'l T *(lf lf 
1 Tx (1 Co - 

-i Tx-t CotCo 
L A l L 2 
■tTxft Co 



lf°lP) 



(29) 



where 





> 4 C , 


a?, 


a?, 


af) 




t <A<(i-2A)A,A(i + A)<i 




(ar£, 


<*2> 


af, 


af) 




2A)A< A<(i-2A)A 


«(1),«(2),«(3),"(4)) = < 


ip-A i 


OL 2i 


a^: 


af) 




2A)A< A< (2A-i)(i-A) 








af- 


af) 


^(2A 


- i)(i-A)<A<(2A-i)(i-A) 








• 


af) 


r(2A 


- i)(i - A) < A < A max , A (2 - A) < i 



(30) 



■min — — (1 — A)(l ~~ A)> A max - 

All possible switching orders of 



^1 — A) A' ® denotes XOR operation. 



according to A are listed in (l30l) and the first and 



the last orders are impossible unless an additional condition is satisfied to make the regions of A valid, 
i.e. A(l + A) < 1 an d A (2 - A) < 1 respectively. Since lf° and lf° are correlated when l Rx = 
(i.e. p 12 / or A / 0), not only lf° and lf° alone but also the identity of lf° and lf° (i.e. lf° 



iCo 
L 2 



or lf°lf° © lf°lf°) can provide information about 1^ and thus can be used to determine CR link 
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availability. This is actually similar to covariance-based detection. For example, if A > (2/3 2 — 1)(1 — Pi) > 

0, lf° and lf° are probably identical when spectrum is unavailable at CR-Rx, i.e. 1 RX = 0, and lf° in 
(|26l ) is then replaced by if ° <g> Furthermore, when A increases and is greater than (2 Pi — 1)(1 — p 2 ), 
the roles of if ° and if ® if ° switch because the identity of if ° and if ° can provide more information 
about 1^ than if ° alone. Alternatively, when A < 0, the results can be similarly explained. In addition, 
it is interesting to note that even if lf° is independent to l Rx (i.e. /3 2 = 1/2), lf° may become helpful 
due to the correlation between lf° and lf°. In the next section, we will further investigate impacts of 
correlation between lf° and lf° on network operation. 

D. Multiple Cooperative Nodes with Limited Statistical Information 

In CRN or self-organizing networks, due to lacking of centralized coordination, each node in CRN can 
only sense and exchange local information. In addition, dynamic wireless channels and mobility of nodes 
make the situation severer and one node can only acquire information within limited sensing duration. 
We can either design systems under simplified assumptions, which may result in severe performance 
degradation, or apply advanced signal processing techniques based on minimax criterion ll42ll . robust to 
outliers in networks, as we are going to do hereafter. 

To derive the optimum Bayesian detection in Proposition 31 we have to acquire joint pmf of spectrum 
availability at K cooperative nodes, which may require long observation interval to achieve acceptable 
estimation error. If we only have up to the /cth order marginal pmf (related to capability of observation), 

1. e Q^, s — 0, 1 according to Lemma H our design criterion becomes minimax criterion, that is, we 

' (s) 

find the least-favorable joint pmf P K , s = 0, 1 such that maximizes Bayesian risk and then conduct the 
optimum Bayesian detection under that joint probability. Therefore, the problem can be formulated as 
follows. 

Proposition 8 (Robust Cooperative Sensing). Cooperative spectrum sensing with limited statistical infor- 
mation Qk K 7 s = 0, 1 becomes ri2?l ) with (P^ \ Pj£ ) replaced by (P^, P^), where 

(PS. P£) = arg max R (p« P<») = arg min ||«,(1 - a )Pg J - aPgUi 

■pi 1 ) -pi 1 -') \ / -p(i) ■pi'-') 

K ' K K ' K 

„ t p( x ) _ (c {1) i-VnW r (1) p (1) ) 

sJ - ^k,K-{^k,K) \^&k,K ~ —k,K—k,K) (31) 

p(0) _ / r (0) \-l/ n (0) p(0) p(0) , 
E-k,K — V^-k,K) \Hk,K ^^K^^K) 

The last equality in the objective function is based on the fact that min(x, y) = (x + y — \x — y\)/2 
and that the sum of probability distribution is equal to one. The result is reasonable because in order to 
minimize the objective function, the likelihood ratio P^ [i]/Pj£ [i] approaches to the optimum threshold 
w(l — a) /a, which induces poor performance of the detector and therefore increases Bayesian risk. 
Furthermore, we could apply Lemma |6] to set the constraints on joint pmf. Since vector norm is a convex 
function, the problem can be solved by well-developed algorithms in convex optimization [1411 . 

In last part, we proposed a simple methodology to select cooperative nodes based on reliability under 
assumption of independent observations. However, in practice, there exists correlation among spectrum 
availability at cooperative nodes and spectrum sensing may change as we showed in Proposition 13 In 
addition, since the statistical information is limited within reasonable observation interval, CR-Tx can 
select cooperative nodes to minimize maximum Bayesian risk by minimax criterion. 

V. Application to Realistic Operation of CRN 

In preceding sections, we only considered single CR link in CRN. However, CRN is not just a link level 
technology if we want to successfully route packets from source to destination through CRs and PS. In 
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the following, we suggest a simple physical layer model for CRN and investigate the impacts of spectrum 
sensing on network operation and the role of a cooperative node playing in CRN, which is impossible to 
be revealed from traditional treatment of spectrum sensing. Since spectrum sensing may not be ideal and 
there exists hidden terminal problem, we further define the true state for PS. 

Definition 6. The true state for PS can be represented by the indicator 




1, PS either does not exist or is inactive 
0, PS exists and is active 



Therefore, with the definition of a, we have 

EUPr(l TiE = 1, 1*" = 1|1 P5 = s) ¥t{1 ps = s) 

ELo Pr (l Tx = l\l PS = s) Pr(l PS = s) 

To connect relations between indicator functions of link availability and realistic operation of CRN, we 
propose a simple received power model. 



A. Received Power Model 

We model the received power from PS and background noise as log-normal distribution, or 10logi (Ps) ~ 
N(/j,s,(Js) and 10log w (P]y) ~ -/V(/z O ;0o)> where erf and a$ are used to quantify the measurement uncer- 
tainty of the received power from PS and noise respectively. In addition, p, is a constant whereas us should 
be varied according to path loss and shadowing. More specifically, let fi S = K — 10alogi (dcR) — ben, 
where K is a constant, d C R denotes distance from CR (either CR-Tx or CR-Rx) to PS as in Fig [2l a 
means path loss exponent, and ben represents shadowing effect. When 1 PS = 1, the received signal only 
comes from noise. However, when 1 PS = 0, the received signal is the superposition of signal from PS and 
noise, which results in addition of two log-normal random variables. We could simply model the received 
power as another log-normal random variable with parameters /icr and <Jcr, and under assumption of 
c*s > we nav e 

{/^o, if Vs<Va- &s 

Us, if Ms > /x + c\s (33) 

(ps + Mo + 0"s)/2, otherwise 

{o"o, if Ms < Mo - crs 

erf, if Ms > Mo + 2a 5 (34) 

^l^rMMs - Mo) + ' 7s+ 3 2 ' 70 , otherwise 

By simulation, the distribution of the simplified model, although not exactly identical to, is close to the 
simulated distribution, especially in terms of mean and variance. It justifies our simplified model. 



B. Spectrum Sensing at CR-Tx and Reception at CR-Rx 

Recall the conditions that CRs can successfully communicate over a link. Assume CR-Tx adopts an 
energy detector in the hypothesis testing (OQ) and there is no interference from co-existing systems. The 
detector can be represented as 

Tx 

Ptx f r Tx (in dB) 

t Tx =0 

where P Tx denotes the received power at CR-Tx and r Tx is a fixed threshold since the detector is designed 
under a given SINR. On the other hand, to successfully receive packets, the SESfR at CR-Rx should be 
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greater than minimum value r] outage as shown in ©. Similarly, spectrum availability at CR-Rx can be 
represented as 



Rx > 



r Rx (in dB) 



where P Rx denotes the received power from PS and noise at CR-Rx. Different from CR-Tx, trx is 
varied according to the received power from CR-Tx. For simplicity, we only consider propagation loss in 
modeling the received power from CR-TX and have r Rx = L — 10alog w (r Rx ), where L is a constant 
and r Rx denotes distance between CR-Tx and CR-Rx. 

We suppose that the measurement uncertainties and hence the received power from PS and noise at CR- 
Tx and CR-Rx are independent. However, to model spatial behavior for CR-Tx and CR-Rx, we consider 
the relation of shadowing between CR-Tx and CR-Rx. Intuitively, the relation should depend on locations 
of CR-Tx, CR-Rx, PS, along with the obstacle size and we proceed based on a linear model 

J max {b Tx {l - 0} , if r Rx cos(6 Rx ) < d Tx , ac . 
0r x — s / n \ j w- 5 ,) 

[0, if r Rx cos(9 Rx ) > d Tx 

where k denotes parameter of obstacle size and 9 Rx is the angle between line segments with starting point 
at CR-Tx and end points at CR-Rx and PS, as shown in Fig. [2] for illustration. In this model, shadowing 
at CR-Rx b Rx linearly decreases with respect to the distance between CR-Tx and CR-Rx r Rx with rate 
inverse proportional to the obstacle size k and is equal to zero when CR-Rx is far apart from CR-Tx 
or PS is located in the middle of CR-Tx and CR-Rx. Additionally, since shadowing parameter achieves 
maximum at CR-Tx, this results in the worst case scenario in spectrum sensing. Finally, from log-normal 
fading distribution, 



a: 



Pr(l PS = 1)Q (fis^a*) Q (^) + Pr(l PS = 0)Q Q ( 



<?Rx, 



Pr(l P5 = 1)Q (^^) + Pr(l PS = 0)Q ( 



where Q(x) denotes the right-tail probability of a Gaussian random variable with zero mean and unit 
variance. 



C. Cooperative Spectrum Sensing 

Under above proposed signal model, the analysis can be easily extended to cooperative spectrum sensing. 
Considering a cooperative node conducting an energy detector, we have 

lC°=l 

Pco i T Co (indB) 

lC"=o 

Furthermore, the correlation due to geography is established similar to (|35T ). We could therefore calculate 
(3 and 7 similar to a, as 

El=o FT ( lTx = h lRx = h lC ° = l\l PS = s) Pr(l P5 = s) 

Ei = o Pr ( lTa; = h lRx = l\l PS = s) Pr(l PS = s) 
ELo Pr (l Ta; = 1, l te = 0, l Co = 0\1 PS = s) Pr(l P5 = s) 
7 ~ ^^ =0 Pr(l r - = = 0|1 P5 = s) Pr(l PS = s) 

With the relation between statistical information {a, (3, 7} and received power model, we can mathemat- 
ically determine allowable transmission region of CR-Tx. 
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D. Neighborhood of CR-Tx 

In Section III and IV, we have developed spectrum sensing under different assumptions and note that 
spectrum sensing depends on the value of a, i.e., spatial behavior of CR-Tx and CR-Rx. Especially, there 
is even a region of a that CR-Tx is prohibited from forwarding packets to CR-Rx and the link from 
CR-Tx to CR-Rx is disconnected (i.e. l hnk = 0). This undesirable phenomenon alters CRN topology 
and heavily affects network performance, such as throughput of CRN, etc. Therefore, we would like to 
theoretically study link properties in CRN and first define the regions of a as follows. 

Definition 7. The set {a\ P T (l link = 0) = 1} is called prohibitive region while {a\ Pr{l link = 1)^0} 
is called admissive region. The boundary between these two sets is called critical boundary of a and is 

denoted by a c . Therefore, {a\l link = 0} = {a|0 < a < a c } and {a\V ink = 1} = {a\a c < a < 1}. 

If a lies in the prohibitive region, the link from CR-Tx to CR-Rx is disconnected. The property and 
the engineering meaning of ac are addressed as follows. 

Lemma 8. «c is a decreasing function with respect to number of cooperative nodes. 

Proof: It is easy to show that for fixed w, a decreases as the threshold of the likelihood ratio test 
w(l — a) /a increases. Therefore, ac can be determine by the largest likelihood ratio. Assume the largest 
likelihood ratio with k — 1 cooperative nodes occurs at i max , imax = ar g max i{Pfe 1 \[^/Pi°-i[^}- 
When the /cth cooperative node enters, let 

Pk — P^lfe* — ° — A^ l)fe _ 1 [l, imax], ■ ■ ■ , = ^-m,k-l[^ ~ 1> ^max], 1^ = 1, = 1) 

Ik — Prllfc 70 — 0|lf° = A^ n jfc _ 1 [l, i m ax\i ■ ■ ■ i lfe-i — A-l n ^ k _ 1 [k — 1, i max \, l Rx = 0, l Tx = 1) 
Then, there are two likelihood ratio with k cooperative nodes, say ith and jth, becoming 

Pj 1} [*] _ p£i[w]& 
Pi 0) M Pi°i[w](l-7*) 

pj^bl _ p£i[w](i-/3 fc ) 

Pfbl P?i[w]7* 

Since either p k + % > 1 or p k + % < 1, one of the ith and the jth likelihood ratio is not less than 
Pfc-it^/Pfc-iW' wn i cn results in lower « c . ■ 
Lemma 9. The following two statements are equivalent: 

1) a > ac 

2) p r (li<nfc = = i) > w /(w + i) 
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Proof: Since a > a c if and only if p r (L mfc = 1) ^ 0, we have 
Pr(l' i,lfc = x|l'* nfc = l) 
= Pr(l Tx = 1, l te = l|l Tx = 1, i te = 1) 
Pr(l fe = 1, i fe = l|l Tx = 1) 

~ Pi(i Rx = i\i T * = l) 

Pr(l^ a: = ljl^ 21 = 1) Pr(i^ a; = = 1 = ]1 

" j2l =Q Pr{l Rx = s\l Tx = 1) Pr(l te = l|l fe = s, 1 T * = 1) 

a Pr(j fl:c =l|l- R3: =l,l Tj: =l) 
l-o Pr(i fl:c =l|l iia: =0,l Ta: =l) 



> 



.. g Pr(l fl3: =l|l- R3: =1,1 T3: =1) 

1-Q Pr(i« a: =l|l Jia: =0,l T:,: =l) 
g ui(l-g) 
1-a a w 



1 + _g^ ™(i-<*) ty + 1 

1— a a 

The inequality holds because the likelihood ratio is greater than w(l — a) /a if l Rx = 1 and x/(c + x) is 
a increasing function with respect to x. Reversely, the conditional probability p r (l /mfe = i|i /mfc = i) [ s 
well-defined if and only if p r (r mfc = 1) ^ 0, which implies a > ac- ■ 
In Lemma |H p r (l hnfe = i|i' mfc = i) could be interpreted as the probability of successful transmission 
in CRN and the weighting factor in Bayesian risk © can be determined by the constraint on the outage 
probability P out = p r (l tofc = 0\l link = 1). That is, if a CRN maintains P out < (, w = (1 - C)/C 
Therefore, the condition that allows CR-TX forwarding packets to CR-Rx (i.e. a belongs to admissive 
region) guarantees the outage probability of CR link. Further considering the proposed physical layer 
models, we can establish and define a geographic region, where CR-Tx is allowed forwarding packets to 
CR-Rx as long as CR-Rx lies in the region. 

Definition 8. Neighborhood of CR-Tx Af is {(tr x , 9r x ) \a > ac} or equivalently becomes {(tr x , 9r x ) | Pr(l hnk 
l\l lmk — 1) > w/(w + 1)}. Coverage of CR-Tx is neighborhood of CR-Tx without PS. 

Please not that the coverage of CR-Tx is defined without the existence of PS and the neighborhood is the 
effective area in real operation coexisting with PS. When we consider a path loss model between CR-TX 
and CR-Rx, coverage becomes a circularly shaped region. However, due to hidden terminal problem as in 
Fig. [Hand [21 where PS is either apart from CR-Tx or is blocked by obstacles, the probability of collision 
at CR-Rx could increase and CR-Tx may be prohibited from forwarding packets to CR-Rx. Therefore, 
neighborhood of CR-Tx shrinks from its coverage and is no longer circular shape. In addition, hidden 
terminal problem is location dependent, that is, PS is hidden to CR-Tx but not to CR-Rx in Fig. \T\ and 
[2l Thus, CR-Rx is possibly allowed forwarding packets to CR-Tx. From such observations, CR links are 
directional and can be mathematically characterized as follows. 

Definition 9. CRi is said to be connective to CRj if CRj is located in the neighborhood of CR it which 



is denoted by lfj lk = 1. Otherwise, l- infc = if CRi is not connective to CRj. 



3 y 

L ij — *- yy^, j-jj — u y CRi iS nOt wim^mv*, iv 

According to above arguments, it is possible that CR-Rx is connective to CR-Tx but the reserve is not 
true. Mathematical conclusion is developed in the following, and is numerically verified in Fig. [6^nd [7] 
in Section VI. 

Proposition 9. Connective relation is asymmetric, that is, for two cognitive radios, CRi is connective to 
CRj does not imply CRj is connective to CRi, or mathematically, l l ^ nk = 1 =£* l l ™ k = 1. 

Proof: We analytically illustrate using Fig. [TJ where CRi lies in the middle of CRj and PS-Tx and 
PS-Tx is hidden to CRj but not to CRi. Let w = 9 to guarantee the outage probability of CR link less 
than 0.1 and let Pr(l P5 = 1) = 0.7, i.e., the spectrum utility of PS is only 30%. If CRi wants to forward 
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packets to CRj (i.e. CRi is CR-Tx and CRj is CR-Rx), CRi can successfully detect the activity of PS 
and Pr(l Tx = 1|1 P5 = 1) » 1 and Pr(l Tx = 1|1 P5 = 0) w 0. Therefore, Ci^ forwards packets to OR, 
only when 1 P5 = 1. In addition, since CRi is located in the transmission range of CRj, CRj is located 
in the transmission range of CRi in a pure path loss model and p r (l i?a ' = l|l Tz = 1, 1 P5 = 1) « 1. 
Applying (l32l) and ©, we have a ~ 1 and l'* nfc = 1. On the other hand, when CRj wants to forward 
packets to CRi, CRj becomes CR-Tx and CRi becomes CR-Rx. Since PS-Tx is hidden to CRj, at CRj, 
the received signal power from PS is below noise power, and [Icr = A*o an d o 2 CR = a\ in (1331 and (|34l ). 
Therefore, Pr(l Tx = 1|1 P5 = 1) « 1 and Pr(l Tx = 1|1 P5 = 0) « 1. That is, CR, always feels the 
spectrum available and intends to forward packets to CR4. However, when 1 PS = 0, collisions occurs at 
CRi and p T {l Rx = l\l Tx = 1, 1 P5 = 0) « 0. Similarly, by ([32]) and ©, we have a « Pr(l P5 = 1) = 

0.7 < to/O + !) = °- 9 and ll j? k = °- ■ 
Proposition [9] mathematically suggests that links in CRN are generally asymmetric and even unidirec- 
tional as the argument in 1146*1 . Therefore, traditional feedback mechanism such as acknowledgement and 
automatic repeat request (ARQ) in data link layer may not be supported in general. This challenge can 
be alleviated via cooperative schemes. Roles of a cooperative node in CR network operation thus include 

1) Extend neighborhood of CR-Tx to its coverage 

2) Ensure bidirectional links in CRN (i.e. enhance probability to maintain bidirectional) 

3) Enable feedback mechanism for the purpose of upper layers 

Since neighborhood increases as a c decreases, by Lemma |«**J the capability of cooperative schemes to 
extend neighborhood increases when number of cooperative nodes increases. Therefore, spectrum sensing 
capability mathematically determine CRN topology. It also suggests the functionality of cooperative nodes 
in topology control [|44l[|45l and network routing lf46l . which is critical in CRN due to asymmetric links 
and heterogeneous network architecture 1146*1 . 

Here, we illustrate impacts of correlation between spectrum availability at cooperative nodes on neigh- 
borhood. Recall Proposition |2l where we considered two cooperative nodes with /3 4 = 7^ p; > 0, i = 0, 1, 
and lf° > From (l30l> . we have 



Therefore, as A increases from A min , ac = ofl increases from and achieves maximum at A = 
(2f3 2 — 1)(1 — Pi). At this point, a c = af = = w(l — 71) /(/Si + w(l — 71)), which is the critical 
boundary with node one alone. If A further increases to A max , ac = decreases to 0. We conclude that 
positive correlation between if and 1^° shrinks the neighborhood, compared to the independent case 
(P12 = 0), unless the correlation is high enough, i.e., A > (2/3 2 — 1)(1 — Pi)/@2 by solving o^|a < 0% |a=o 
according to (l36l) . 

If one CR has larger neighborhood area, it is expected to be connective to more CRs and to have higher 
probability to forward packets successfully and higher throughput of CRN accordingly. The result offers 
a novel dimension to evaluate cooperative nodes. That is, different from criterions in link level, such as 
minimum Bayesian risk or maximum reliability as we mentioned in last section, maximum neighborhood 
area is a novel criterion to select the best cooperative node from the viewpoint of network operation. 

Proposition 10. (Optimum Selection of Cooperative Node) For a CRN with a constraint on the outage 
probability P out < C there are one CR and K cooperative nodes, indexed by k. The best cooperative 
node for the CR under maximum neighborhood area criterion is 



where A/a(&) represents neighborhood area of the CR with the aid of the kth cooperative node. 

In CRN, CRs could act as relay nodes to relay packets to the destination. Assume the destination is in 
the direction 9 of a CR with respect to PS. It is intuitive for the CR to forward packets to the direction 




(36) 



k op t = argmaxA/yi(A:) with w = (1 — C) /C 



k 



(37) 



18 



around 9, say 6±e. Let No± e = {M\9 Rx E (8 — e,6 + e)} and then the best cooperative node may become 
(1371) with J\f replaced by J\f e±e . 

VI. Experiments 

A. General Spectrum Sensing 

1 ) Spectrum Sensing Performance: The performance of spectrum sensing, measured by Bayesian risk 
©, is plotted by Bayesian risk versus the probability of spectrum availability at CR-Rx a in Fig. [3] 
We set the weighting factor w = 9 (w is defined in ©) to guarantee the outage probability of CR link 
less than 0.1. Larger Bayerian risk represents worse performance because spectrum sensing induces more 
possibility of collisions at CR-Rx or of losing opportunity to utilize spectrum. We see that traditional 
spectrum sensing without considering spectrum availability at CR-Rx (i.e. l Rx ) has large Bayesian risk 
when a becomes small because collisions usually occur when CR-Tx determines link availability only by 
localized spectrum availability at CR-Tx (i.e. l Tx ). On the other hand, by considering l Rx in our general 
spectrum sensing © with known a, Bayesian risk decreases when a is less than 0.9, which is the critical 
boundary of a, etc, and CR-Tx is prohibited from forwarding packets to CR-Rx. Therefore, risk occurs 
due to loss of opportunity to utilize spectrum. 

However, in practice, a is unknown and needs to be estimated by Lemma HI We set observation depth 
(i.e. duration) L = 15 and show expected Bayesian risk of inference-based spectrum sensing © with 
respect to observed sequence l Rx [n — 1], ... , \ Rx [n — L\. The performance degrades around a c and even 
worse than that of traditional spectrum sensing when a > a c because the estimation error may cause the 
estimated a © to across ac and results in different sensing rules; however, it is close to the performance 
with known a. This verifies the effectiveness of inference-based spectrum sensing. 

Fig. |3] also shows Bayesian risk of cooperative sensing © under different sensing capability of a coop- 
erative node, i.e. reliability M R . We assume statistical information {a,P,j} can be perfectly estimated. 
The performance curve is composed of three line segments as in © and shows performance improvement 
in the middle segment due to the aid of the cooperative node. However, in the right and the left segments, 
the cooperative node becomes useless and the performance is equal to that of non-cooperative sensing. 
In comparison of sensing capability of a cooperative node, the one with larger reliability is expected to 
have higher correlation between spectrum availability at CR-Rx and the cooperative node (i.e. l Co ) and to 
provide more information about l Rx ; therefore, it achieves lower Bayesian risk and lower a c (i.e. larger 
admissive region). In addition, when /3 + 7 = 1 (/3 = 0.8, 7 = 0.2), l Rx and l Co are independent and the 
performance is identical to that of spectrum sensing with known a in Fig. |3] 

2) Impacts of Correlation between Spectrum Availability at Cooperative Nodes: We next investigate 
performance of spectrum sensing with two cooperative nodes with respect to the correlation between 
spectrum availability at these two nodes (i.e. lf° and pi 2 . We set f3\ = 71 = 0.75, /3 2 = 72 = 0.7 
as the scenario in Proposition |7] and depict Bayesian risk in Fig. HI Generally speaking, Bayesian risk 
decreases (increases) as p 12 increases in high (low) a. We also observe that ac decreases when number 
of cooperative nodes increases and a c increases when p 12 increases unless the correlation is high enough. 
For example, if two cooperative nodes are close in location and p i2 = 0.8, ac is less than that when 
p\2 = 0. It is also interesting to note that there exists a region of a such that the identity of lf° and 
I2 instead of if determines CR link availability as in (|29| ) and Bayesian risk is less than that with if 
alone. 

The results further suggest trade-off between performance in link layer and network layer when we 
select cooperative nodes. That is, for one CR link with a > 0.5 (e.g. CR-Rx is close to CR-Tx or spectrum 
utilization of PS is low), large p\% is preferred to achieve low risk. However, from network perspective, 
to achieve high number of CR links that are admissive to CR-Tx (i.e. to achieve large neighborhood and 
low ac) and thus high throughput of CRN, small p 12 is preferred. 
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3) Robust Spectrum Sensing: For multiple cooperative nodes, with six nodes in our simulation, we 

(s) 

show the performance of cooperative spectrum sensing with limited statistical information Q k K , s = 0,1 

due to limited sensing duration. We first find least-favorable joint pmf Pw, s = 0, 1 by (f3TT ) and then 
compute corresponding Bayesian risk, which is shown in Fig. [5] under different order of known marginal 
pmf k (i.e. capability of observation). That is, CR-Tx only acquires pmf of k out of six cooperative nodes. 
The risk is compared to that with the optimum sensing rule (T23T) and that with assumption of independent 
spectrum availability ( 128~T ). Obviously, Bayesian risk decreases and approaches to that in the optimum 
case as the order of known marginal pmf k increases because more information is acquired to generate 
Pjw, s — 0, 1 closer to the true one , s — 0, 1. We observe that when the order k is greater than 3, 
robust spectrum sensing outperforms the case of traditional independence assumption. Therefore, if CR-Tx 
would like to select six cooperative nodes, CR-Tx only requires statistical information about spectrum 
availability among three out of six cooperative nodes, i.e. Qjjg, s = 0, 1, to achieve better performance 
than the case according to reliability criterion. 

B. Neighborhood of CR-Tx 

1) Without Obstacles: In Fig. |6j we illustrate neighborhood of CR-Tx ("+" in the figure) without 
blocking. The neighborhood boundary with/without a cooperative node ("o" in the figure) is depicted by 
a thick and a thin line respectively. The parameters are set as follows: fi = 0, <Tq = 1, a 2 s = 8, K = 10, 
a = 3, L = 3, t Tx = T Co = 3, and Pr(l P5 = 1) = 0.6. 

In Fig. Oa), PS ("*" in the figure) is placed near to CR-Tx ((0.7, 0)). We observe that CR-Tx almost 
perfectly detects the state of PS, i.e., Py{1 Tx = 1\1 PS = 1) « 1 and Pr(l Tx = 1|1 P5 = 0) « 0, and 
a p r (l Rx = l|l Tx = 1, l ps = 1) by (|32l) . Therefore, neighborhood of CR-Tx approaches to its 
coverage and the cooperative node is not necessary in this case from a viewpoint of network operation. 
However, when PS is apart from CR-Tx ((1.7,0)) as we shown in Fig. [T] the neighborhood at PS side 
shrinks and is no longer circularly shaped because PS is hidden to CR-Tx and hence probability of collision 
at CR-Rx increases when 1 P5 = 0. Fig. |6fb)~(d) illustrate the neighborhood under different locations 
of the cooperative node. We observe that neighborhood area decreases when the cooperative node moves 
away from PS and there even exists a region where cooperative sensing can not help. Therefore, the 
cooperative node in Fig. Ob) is the best among these three nodes according to maximum neighborhood 
area criterion in Proposition [10l 

We present an example of existence of unidirectional link in CRN. In Fig. |6fb), assume one CR is 
located at (1,0). Obviously, the CR-Tx is not connective to the CR and therefore is prohibited from 
forwarding packets to the CR. However, by Fig. [6ta), the CR is connective to CR-Tx, which makes the 
link unidirectional (only from the CR to CR-Tx). As Proposition |9l this also shows asymmetric connective 
relation even under rather ideal radio propagation. With the aid of a cooperative node located at (0.4, 0.3), 
the link returns to a bidirectional link. 

2) With Obstacles: Alternatively, we consider effects of shadowing due to blocking, as we illustrated 
in Fig. [2l We set shadowing parameters bx x = 25 and parameter of obstacle size k = 0.3 and 0.7 in 
Fig. I2ta)(b) and Fig. [7jc)(d) respectively. We observe that small obstacles size (i.e. small k) can result in 
more substantial shrink of the neighborhood, compared to large obstacles size (i.e. large k). The reason 
is: if k is small, only a small region around CR-Tx falls in deep shadowing and the state of PS can be 
successfully detected outside that region. Therefore, this leads to high probability of collision at CR-Rx 
as 1 P5 = 0. On the other hand, if k is large, CRs are likely separated from PS by obstacles, which results 
in large "distance" between CR and PS. Here, "distance" is measured by received signal power [fT6ll[[36ll . 
In comparison of the capability of a cooperative node, the one in small k has good capability to recover 
the neighborhood to its coverage even when the node is at opposite side of PS. However, for large k, 
the cooperative node may also be in deep shadowing and becomes useless to recover neighborhood of 
CR-Tx. 
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VII. Conclusion 

In this paper, we showed that CR link availability should be determined by spectrum availability at both 
CR-Tx and CR-Rx, which may not be identical due to hidden terminal problem (Fig. __ and [2]). In order 
to fundamentally explore the spectrum sensing at link level and its impacts on network operation, we 
established an indicator model of CR link availability and applied statistical inference to predict/estimate 
unknown spectrum availability at CR-Tx due to no centralized coordinator nor information exchange 
between CR-Tx and CR-Rx in advance. We therefore expressed conditions for CR-Tx to forward packets to 
CR-Rx under guaranteed outage probability. These conditions, along with physical channel models, define 
neighborhood of CR-Tx, which is no longer circularly shaped as coverage. This results in asymmetric or 
even unidirectional links in CRN, as we illustrated in Section VI. The impairment of CR links can be 
alleviated via cooperative scheme. Therefore, spectrum sensing capability determines network topology 
and thus throughput of CRN. Several factors with impacts on spectrum sensing are analyzed, including: 

1) Correlation of spectrum availability at cooperative nodes 

2) Capability of observation at CR-Tx (i.e. available statistical information at CR-Rx) 

3) Locations of cooperative nodes and environment (i.e. obstacles) 

Furthermore, limits of cooperative scheme were also addressed in link level and network level. In 
addition, to measure sensing capability and then to select cooperative nodes is an important issue because 
we would like to minimize information exchange to increase spectrum utilization. Criterions from link 
level (maximum reliability or minimum Bayesian risk) and network level (maximum neighborhood area) 
perspectives were accordingly proposed. We numerically demonstrated existence of trade-off in designing 
systems in different layers. In addition, robust spectrum sensing was proposed to deal with local and 
partial information due to no centralized coordination and limited sensing duration in CRN. More useful 
results in CRN extended from this research can be expected in future works. 
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Fig. 1. Hidden terminal problem. CR-Rx lies in the middle of CR-Tx and PS-Tx and PS-Tx is hidden to CR-Tx. 
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Fig. 2. Activities of PS-Tx is blocked by obstacles to CR-Tx but not to CR-Rx. 
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Spectrum Sensing Performance, w=9 




Fig. 3. Performance (Bayesian risk) of spectrum sensing schemes with respect to the probability of spectrum availability at CR-Rx a. The 
critical boundary of a, ac, is represented by • under different sensing schemes . 
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Fig. 4. Performance of spectrum sensing with one/two cooperative node(s) with respect to the probability of spectrum availability at CR-Rx 
a. The transition and the value of the critical boundary of a, ac, under different correlation of Node 1 and Node 2 pi2 is also shown. 
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Order of Known Marginal pmf 



Fig. 5. Performance of robust spectrum sensing under different order of known marginal pmf k, that is, pmf of k out of six cooperative 
nodes is available at CR-Tx. The performance is compared to the spectrum sensing with full statistical information and with traditional 
independence assumption. 



25 




Fig. 6. Neighborhood of CR-Tx without obstacles. CR-Tx can only be allowed forwarding to CR-Rx located within the bounded region. 
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Neighborhood of CR-Tx, With Obstacles, w=4, b =25 
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Fig. 7. Neighborhood of CR-Tx with obstacles. Effects of small (n = 0.3 in (a)(b)) and large (k = 0.7 in (c)(d)) obstacles are compared. 



